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The smallness of the 1 — 3 lepton mixing angle #13 and of the neutrino mass-squared-difference 
ratio Am? 2 /Airi23 can be understood as the departure from a common limit where they both vanish. 
We discuss in general the conditions for realizing the mass degeneracy of a pair of neutrinos and 
show that the vanishing of a CP violating phase is needed. We find that the discrete quaternion 
group Q of eight elements is the simplest family symmetry which correlates the smallness of Ami 2 
to the value of #13. In such a model we predict 0.12 < sin 613 < 0.2 if the ordering of the neutrino 
mass spectrum is normal, and sin #13 < 0.12 if it is inverted. 
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Introduction. Contrary to quarks and charged leptons, 
the three neutrinos are known to have a mass spectrum 
with a weak hierarchy, possibly quasi-degenerate. In par- 
ticular, the two mass eigenstates v\ and vi in the "solar 
pair" are very close in mass in the case of inverted order- 
ing and also in the case of normal ordering, as long as 
the yet unknown absolute mass scale is larger than ~ 0.02 
eV. Indeed, a global fit of neutrino oscillation data gives 

Am 2 2 = m\ - m\ = (7.9±J;g) ■ 10~ 5 eV 2 , 

Am^ 3 ee ml - m\ = ±(2.6 ± 0.6) • 10" 3 eV 2 , 

sin 2 9 12 = 0.30±g-S° , (1) 

sin 2 9 23 = 0.50t° H , 

sin 2 6»i3 < 0.040 , 

where we took the best fit values and the 3cr intervals 
from the last update of Ref.pQ (v5). 

The smallness of the 1 — 2 mass splitting compared 
to the 2 — 3 "atmospheric" splitting may be explained 
as the departure from a symmetric limit where 1—2 
mass degeneracy holds. Other small parameters of the 
lepton flavor sector may also be interpreted as deviations 
from the same limit, such as the 1 — 3 mixing angle, the 
deviation of the 1 — 2 and possibly 2 — 3 angles from 
the maximal value 7r/4, and the mass ratios {m e /m^, 
m^/m T ) of charged leptons. 

In this paper we analyze the neutrino mass structures 
corresponding to a quasi-degenerate pair of states and 
the phenomenological correlations with other lepton fla- 
vor parameters, in particular a non- vanishing 1—3 mix- 
ing, whose measure is the objective of an extensive ex- 
perimental program [2J. We also search for the simplest 
flavor symmetries which can be used to realize such a 
mass degeneracy. The discrete quaternion group of eight 
elements Q is identified as the most suitable for this pur- 
pose and a complete model is constructed. 

The group Q (sometimes called Q s or Q4) was 
introduced in [3] to build a model of quark and lepton 
masses and mixing. Discrete subgroups of quaternions 
with unit norm, i.e. SU(2), were already used in [4] 
to suppress the neutrino mass while allowing for large 



neutrino magnetic moments. They were also discussed 
as flavor symmetries in a series of papers by Frampton 
and collaborators [5j. A specific model based on the 
quaternion group of 12 elements Qe was proposed as well 
[6]. Models were built [7] using the binary tetrahedral 
group T' (double covering of the tetrahedral group A±), 
which is also a discrete subgroup of quaternions with 
24 elements. The group T' was recently employed to 
accommodate tri-bi-maximal mixing [5]. Note that A4 
is not a subgroup of T' , but Q is. Geometrically, SU(2) 
is isomorphic to the hypersphere in four dimensions, the 
8 elements of Q form the 8 vertices of the perfect hy- 
peroctahedron (dual of the hypercube), whereas the 24 
elements of T" form the 24 vertices of the hyperdiamond 
(which is self-dual). 

Mass matrix of 2 degenerate neutrinos. Let us 

begin considering the Majorana mass matrix m v for two 
neutrino states. Under the requirement of mass degener- 
acy (mi = m2 = m), m v can be written in full generality 
as 

„ „ 1 a b \ ( cos 29 sin 29 \ 
D ^ D * = { b -a ) = { sin20 -cos20 ) m ' ® 

where D$ — diag(e 1 ^' 1 ,e^ 2 ) is a diagonal matrix of 
phases, a and b are real and positive, m = \J d 1 + b 2 
and tan 29 = b/a. In terms of the m v matrix elements, 
the requirement of mass degeneracy is equivalent to 2 
conditions, |m n | = \m 2 2\ and arg(m 11 m 2 2/?Ti 1 2) = tt- 

The simplest cases b = or a = have often been 
discussed, while the general case was studied in just a 
few interesting papers, as a prototypical example of a 
pseudo-Dirac neutrino mass matrix [9j I10| 111] . Here we 
analyze in detail how this matrix structure relates to the 
physical observables, with special attention to the effects 
of possible CP violating phases. 

The diagonalization of the matrix in eq.Q presents 
some subtleties, which turn out to be important to un- 
derstand the effect of small perturbations responsible for 
mi =/= mi. Notice first that the two neutrino masses mi,2, 



2 



and the moduli of the mass matrix elements, a and b, 
are physically well-defined quantities (they can be mea- 
sured, at least in principle). This determines uniquely 
the parameter 9, which one is tempted to identify with a 
physical mixing angle 9i2 between v\ and v^. However, 
in the degenerate limit there is no mixing angle responsi- 
ble for neutrino oscillations: m v m\ — diag(m 2 ,m 2 ). To 
settle this apparent contradiction let us rewrite m„ in the 
standard parameterization, 



D' 



C12 
-S12 



S12 
C12 



me 2lp 




C12 
S12 



-812 
C12 



D' 



(3) 

where D'^ is a diagonal matrix of phases, C12 = cos #12, 
s\2 = sin 0i2 and p is the relative Majorana phase be- 
tween the two mass eigenstates, varying between and tt. 
With some easy algebra one finds sin 29 = sin psin 2(?i2- 
Only this combination has physical meaning (it is mea- 
surable), while the mixing angle #12 and the Majorana 
phase p cannot be determined uniquely (a similar discus- 
sion can be found in |12j). 

The important consequence is that different small 
perturbations which generate Am\ 2 may select very 
different values for #12, which is the crucial parameter 
to determine the oscillation probability. Let us consider 
a positive mass e <C a,b. If m v diagonal entries are 
corrected as a — > (a — e) and —a — * (—a — e), then 
Am 2 2 = 4eV a 2 + b 2 but the mixing angle is a free 
parameter, tan20 12 = b/a, and the two neutrinos have 
opposite CP-parity, p = tt/2. If instead m v off-diagonal 
entries are corrected as b — > (b + ie) , then Am\ 2 = 4ea 
and the mixing angle is maximal, 9i2 = tt/4, but the 
Majorana phase is a free parameter, sin p s» b/^/a 2 + b 2 . 

Quaternion model. We now search for family symme- 
tries that can lead to the structure in eq.Q. The equal- 
ity I mil I = 1^22 1 cannot be explained by an Abelian 
symmetry, since in this case each lepton family would 
transform independently under the action of the symme- 
try group, so that equalities among independent mass 
matrix elements cannot be justified. Hence the two 
quasi-degenerate neutrino families should sit in a two- 
dimensional irreducible representation (2-dim irrep) of a 
non-Abelian group. Let us assign, therefore, the three 
Standard Model (SM) lepton doublets as follows: 



Li 
L 2 



2 , L 3 



(4) 



where the 1-dim irrep is not necessarily the singlet in- 
variant under the symmetry group. To realize the struc- 
ture in eq.Q, one needs two invariants, a(L\L\ — L2L2) 
and 2&L1L2, and at the same time the combination 
L\L\ + L2L2 should not contribute. We found that this 
2x2 pattern may be obtained by using the 2-dim irrep 
of any of the three smallest non-Abelian groups, 53 [IB] . 
D4 P3] and Q. However, in the S3 case it is not possible 



to maintain this pattern in a complete model with three 
families. In the D4 case the combination + L2L2 is 
a group invariant and cannot be discarded without extra 
assumptions. 

We therefore focus on the smallest quaternion group 
Q. All the details on the group structure, the character 
table and our conventions for the irreps and their tensor 
products can be found in [I3J. For our purposes here it 
is sufficient to recall that Q has four 1-dim irreps, 1 ++ , 
10 , 1 h and 1 , with tensor product rules made ob- 
vious by the superscripts, and one 2-dim irrep 2. The 
product of two Q-doublets (ipi 02) T , (xi X2) T ~ 2 goes 
as follows: 



("01X2 - 02Xi) 
^1X2 + 02Xi) 



(V'iXi - ^2X2) 

("01X1 +^2X2) 



1+ 
1 



(•5) 



The charged lepton masses arise from the Yukawa cou- 



pling yijkLi 



. where d>h = 



,(/> )k are Higgs doublets 



with vacuum expectation values (VEVs) (tpj,) = v%. The 
neutrino masses arise from the Majorana-type Yukawa 
coupling fijkLiLjAk where A k = (A ++ ,A+,A°) k are 
Higgs triplets with (A") = Uk- These VEVs are natu- 
rally small when the triplets are super-heavy, by virtue 
of the type II seesaw mechanism. 

Let us consider the following Q assignments: 



Li 
L 2 



2, L 3 ~l++, e?~l-+, 1— , 

2, 03-1—, Ai~l-+ 1+-, 1++ 



(6) 



Then the charged lepton and neutrino mass matrices have 
the following structure: 



mi 




-J/2«l Z/3«l 
y2"2 V3V2 
V4V 3 

fiui 
-f 2 u 2 
/3M3 



(7) 



The 1 — 2 sector of m u has the desired form of eq.Q. 
Therefore, in the CP conserving case, the 1 — 2 mass de- 
generacy is realized. Moreover, when Q is broken by the 
Higgs doublet VEVs in the direction V\ = 0, the unique 
off-diagonal element in the charged lepton mass matrix 
is (mz)23- Therefore the 1 — 2 mixing comes entirely from 
m u , the 2 — 3 mixing entirely from mi and one predicts 
9 13 = 0. The scalar potential, the generation of VEVs 
and their alignment are discussed in the Appendix. 

A possible extension of the Q symmetry to the quark 
sector is discussed in [3], where the phenomenological 
constraints on the extra Higgs doublets are also esti- 
mated. Notice that, since charged leptons mix only in 
the 2 — 3 sector, there are no flavor-changing neutral cur- 
rents (FCNCs) involving the electron, which are strongly 
constrained experimentally. Several possible tests of FC- 
NCs in the p — t sector are discussed e.g. in [T5] . 
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Before performing a detailed analysis of lepton masses 
and mixing angles associated with the matrices in eq.Q 
and their small perturbations, we would like to stress that 
the same pattern is maintained in many possible variants 
of the Q model, with a different symmetry breaking sec- 
tor, different field assignments and/or a different type of 
seesaw. 

First, one may dislike the presence of multiple Higgs 
doublets at electroweak scale, because of sizable FCNC 
effects, a harsher hierarchy problem (more than one fine- 
tuning), worsened gauge coupling unification, etc. Of 
course all such worries are based on some amount of 
theoretical prejudice. In any case, one can rephrase 
the flavor model above in terms of a unique Higgs 
doublet <f> invariant under the family symmetry, then 
adding SM singlets (pi (flavons) charged under Q as 
(2, 1 ). They enter charged lepton Yukawa couplings 
as yijkLie^((pk/A.)(f> where A > (cpk) is some cutoff scale. 
In this way one can reproduce the same mass matrix 
structure as before, while maintaining the SM particle 
content only at electroweak scale. In this context, super- 
heavy triplets may also be eliminated. One may think of 
neutrino masses originating from the effective operator 
fijkLiLj((pk/ K)(cj)* <j)* / h), with flavons ipk charged under 
Q as 1 h and 1+ . 

Second, a different structure of mi leading to the same 
mixing pattern as in eq.(|7]) can be obtained also when 
charged lepton singlets e\ and lepton doublets Li both 
transform as (2, 1 ++ ), which is required in left-right 
symmetric extensions of the SM. In this case, adding 
04 ~ 1 H , one finds 

/ V3V3 - VaVa -yiv 2 \ 

"ii= y 3 v 3 + y 4 v 4 y 2 v x . (8) 

\ -yiV2 yivi o J 

When Q is broken by Higgs doublet VEVs in the direc- 
tion V2 — 0, one can accommodate m e ^ T and, at the 
same time, large (maximal) 2 — 3 mixing. 

Third, the neutrino mass matrix in eq.(|7]) can be 
derived with just a little bit more effort even in the 
context of type I seesaw, by introducing neutrino 
singlets i/f ~ (2, 1 ). The neutrino Dirac mass 
matrix tub comes from the Yukawa coupling y^L^c/), 

where <j) = (</> + ,0o) is a Higgs doublet transforming as 
1 , so that mj = diag(x,x,y). The neutrino singlet 
mass matrix Mr comes from the coupling fijkvfVjSk, 
where Sk are Higgs singlets with superheavy VEVs, 
transforming as (1 ^ l" 1 , 1 ++ ). In this way one finds 
that Mr, M^ 1 and m„ = — m^-M^m^ all have the 
same structure as m v in eq.([7|. A similar model with 
type I seesaw which realizes the same form of m„ by 
means of a SU(2) x U(l) family symmetry can be found 
in [H]. 

Correlations among observables. The neutrino 
mass matrix in eq.([7]) represents an interesting limit: 
#13 = and either Am] 2 = 0, in the case of no CP viola- 
tion (that is, when (/2U2)/(/iMi) is real), or 6*12 = tt/4, 



when a nontrivial CP violating phase is present. In both 
cases solar oscillation data (see eq.Q) call for a pertur- 
bation to this matrix structure. In the CP-conserving 
case, the perturbation leads in general to a correlation 
between the values of Am^ and #13. In the CP-violating 
case, the correlation will be between the deviation of 1 — 2 
mixing from maximal and nonzero 1 — 3 mixing. Such 
correlations will be probed in future searches of #13 [2]. 

What may be the origin of such perturbations? A first 
possibility is that the structures in eq.([7]) are modified 
by radiative corrections from some large scale down to 
the electroweak scale. This scenario can be justified as- 
suming that the family symmetry is broken at the large 
scale (by the VEVs of flavon fields) and that m v runs 
below that scale as in the SM (or the MSSM). The ra- 
diative corrections to m v for matrix structures similar to 
those considered here have been studied in great detail 
in (T2J [T71 |T5]. It is found that it is possible to gener- 
ate radiatively both Amf 2 and #13, but it is problematic 
to obtain the present values of parameters in the Large 
Mixing Angle MSW region. We will not consider this 
radiative possibility in the following. 

A more straightforward way to introduce a perturba- 
tion is to add extra Higgs multiplets, with different Q as- 
signments, which may provide a sub-dominant contribu- 
tion to m„. First, consider a Higgs triplet A4 ~ 1 . Its 
VEV gives an equal contribution to the 11- and 22-entry 
of m v . In the CP conserving (violating) case, a small per- 
turbation of this type generates Arnf 2 (6> 12 tt/4) 
but does not affect 6*13 = 0. Each observable is repro- 
duced by a different parameter and therefore no corre- 
lations are predicted. In particular one cannot tell this 
scenario from any other model with #13 = 0. 

The most interesting scenario is obtained adding, in- 
stead, (A4, A 5 ) ~ 2, with 114 = (same direction in 
group space as for the Higgs doublets (0 1; </> 2 ) ~ 2 with 
Vi = 0, see the Appendix). Then the neutrino mass ma- 
trix has the form 

fa b d\ 
m u = \ b -a , (9) 
\d c J 

where d is proportional to u§. Let us recall that we are 
working in a basis where mi contains an arbitrary 2 — 3 
mixing 6 l 23 , so the 2 — 3 mixing #23 in m„ is not required 
to match the observed value of 623 — 6 l 23 + 9% 3 . Neutrino 
mass matrices with one zero element and two indepen- 
dent nonzero elements equal to each other are a typical 
outcome of models based on the family symmetry Q [3] . 
All possible matrices with this feature in the basis where 
mi is diagonal were analyzed in |191 . 

CP-conserving case. Let us diagonalize m„ in eq.([9]) 

in the case where all matrix elements are real (without 
loss of generality one can take a, b and d positive). In 
the limit d = e <C a, b, |c|, defining m = yja 2 + b 2 and 
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expanding in e, the three neutrino masses are given by 



mi 



to - 



e 2 (m + a) 
2to(to — c) ' 



m 2 



to 3 



e 2 (c+ a) 



e 2 (a — to) 
2to(to + c) 



(10) 

(strictly speaking these equations hold only as long as 
e <C \m ± c|). The mass squared differences are then 
easily derived: 



Ato 2 2 

^23 = '"-3 

The mixing angles are 
e(c + a) 



2m(a + c) 



„2 



TO 
„2 



(11) 



sin f 13 



C — TO 

sin "" 



23 



tan 2 #12 

r>2 _ ryy-> 2 



(12) 



The almost maximal 2 — 3 mixing arises from the charged 
lepton sector. The ordering of the mass spectrum is nor- 
mal (inverted) for c 2 — to 2 > (< 0). Since tan 2 6 , i2 < 1, 
solar neutrino data require Ato 2 2 > 0, that is a + c > 
(< 0) for the case of normal (inverted) ordering. Finally 
and most importantly, the value of #13 is correlated with 
the other observables: 



sin 2 6 



13 



1 ATO 2 2 

2 Ato^ 3 



^cos 



2(9i 



TO3\ 
TO / 



(13) 



In the case of normal ordering, 777.3 > m and #13 decreases 
by increasing the absolute neutrino mass scale, with a 
lower bound 

sin0i 3 |„ orma i > cos 6»i2 Aml 2 /Amj 3 « 0.15 . (14) 

In the case of inverted ordering 7773 < — tocos2#i2 and 
#13 increases by increasing the absolute mass scale, with 
an upper bound 



Sin C713 I inverted 



< sin6»i2\/A 



IAj 



z 23 I 



0.10 . (15) 



These correlations are the main predictions of our model, 



'12 



and 9 



13- 



which links the smallncss of A7 

We performed a numerical analysis for the mass matrix 
in eq.([9]), scanning over the values of a, b, c and d. The 
prediction for sin #13 is shown in FigjT] as a function of 
the neutrino mass-squared-difference ratio Atoi 2 /Ato 2 3 . 



The bounds (14 1 and (15), which assume best fit val- 



ues of the measured parameters, are slightly relaxed but 
hold qualitatively: we find sin 0i 3 \ norma i > 0.12 and 
sin6*i 3 | int , erte( j < 0.12. 

The prediction for sin #13 is shown in Fig|2] as a func- 
tion of the absolute neutrino mass scale. The lower den- 
sity of allowed points for larger absolute mass scales in- 
dicates that a quasi-degenerate spectrum requires some 



sin 13 



Normal Ordering 




Aitv 
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Am, 



sin © 13 



Inverted Ordering 




Amf 



0.01 0.02 0.03 0.04 0.05 0.06 
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FIG. 1: The 1 — 3 mixing as a function of Am^/Am^ for 
the matrix m v in eq.|9]l, assuming no CP violation. The 
upper (lower) panel corresponds to the case of normal (in- 
verted) ordering of the mass spectrum. We imposed the con- 
straints 0.24 < sin 2 0i2 < 0.40 (3a allowed interval) as well 
as ^ i mi < 1 eV. The dashed rectangle indicates the 3a al- 
lowed region defined by 0.022 < Am^/Am 2 ^ < 0.045 and 
sin 6i3 < 0.2. The gray points are those where d is signifi- 
cantly smaller (less than one third) than a, b and c|. 



fine-tuning in the input parameters a, b , c and d. In 
the normal ordering case the experimental constraint 
sin 6*i3 < 0.2 implies a lower bound mi > 0.02 eV. In 
the inverted ordering case sin #13 vanishes for 7773 ~ 0.02 
eV. A lower bound on sin #13 holds also in the case of 
inverted ordering if 7773 is sufficiently large. This can be 
probed in neutrinoless 2(3 decay searches, which can mea- 
sure the effective mass parameter 777 ee = |mn| = a. In 
the limit of small d one has m ee tocos2$i 2 ~ m/2. 
If m ee ^ 0.04 eV is found, our model predicts a lower 
bound sin #13 > 0.05. 

One should notice that the matrix in eq.([9| may ac- 
commodate data even when d is not much smaller that 
the other parameters, as shown by the black region in 
Figsjl] and |] Indeed, Ato 2 2 and #13 both vanish not 
only in the limit d — > 0, that was studied above, but also 
in the limit c + a — »• 0. When a + c = one has 



ii ; 2 = ±777 = ±Vtt 2 + b 2 + d 2 , 7773 = —a 



tan #23 = — r i tan2#i2 = 



Vd 2 + b 2 



(16) 



Therefore the ordering of the mass spectrum is inverted 
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sin 813 Normal Ordering 
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FIG. 2: The 1 — 3 mixing as a function of the absolute mass 
scale for the matrix m v in eq.Q, assuming no CP violation. 
The upper (lower) panel shows the allowed region in the mi — 
sin 813 (ma — sin #13) plane for the case of normal (inverted) 
ordering of the mass spectrum. We imposed the constraints 
0.022 < Am? 2 /Am^ < 0.045 and 0.24 < sin 2 12 < 0.40 {3a 
allowed intervals). The dashed line indicates the 3a upper 
bound sin S13 < 0.2. The gray points are those where d is 
significantly smaller (less than one third) than a, b and |c|. 



sin 813 Normal Ordering 




sin 8 13 Inverted Ordering 




FIG. 3: The allowed region in the plane sin 2 612 — sin #13 for 
the matrix m„ in eq.|9]l with the most general CP violating 
phases. The upper (lower) panel corresponds to the case of 
normal (inverted) ordering of the mass spectrum. We imposed 
the constraints 0.022 < Am^/Am^ < 0.045 (3<r allowed in- 
terval) as well as m; < 1 eV. The dashed rectangle indi- 
cates the 3a allowed region defined by 0.24 < sin 2 9 12 < 0.4 
and sin #13 < 0.2. The gray points are those where \d\ is 
significantly smaller (less than one third) than aj, \b\ and cj. 



and both large angles can be accommodated (there is no 
need of large 2 — 3 mixing in mi). Defining e = a + c and 
taking the limit |e| <C a, b, d, we find 



tion from maximal 1 — 2 mixing. One should diagonalizc 
M 2 + H 2 + |e| 2 ab* -a*b ae*+c*e 



m u m! 



\b\< 



c*b 



Am 2 12 



2d 2 me 



b 2 



sin 9 



d 2 ' 
bde 



13 



(b 2 + d 2 ) 3 / 2 



1 Am 2 2 sin 26>i 2 

2 Am 23 tan 6*33 



(17) 



'12 



When 6^0, 613 vanishes but at the same time Ar 
is nonzero. This is why in the case of inverted ordering 
there is no lower bound on sin # 13 (see lower panel in 
Figs Jl] and. |2). 

CP-violating case. Let us consider the neutrino mass 
matrix in eq.(|9]) in the general case of complex matrix 
elements. In the limit d = one has 6> 12 = 7r/4, 13 = 0, 
Am 2 2 = 4Im(a6*) and Am| 3 = |c| 2 - \a\ 2 - \b\ 2 . Notice 
that, in the presence of non-trivial phases, the smallness 
of ATO 2 2 /Am 2 3 is accidental. 

When d = e is small, that is, |e| <g; \a\, \b\, |c|, a small 
9%3 is generated and its value is correlated to the devia- 



At leading order, we find 



sin 013 



A j 



1 



\P»(a,b,c)\ 



Am 2 2 



Am 2 3 



(18) 



(19) 



where P 2 is a lengthy expression quadratic in a, 6 and c. 
The form of sin #13 is analog to the one in eq.(12l. The 



deviation from maximal 1 — 2 mixing (sin 2 12 = 1/2) can 
be of order one for |e| 2 > Amf 2 . 

We performed a numerical analysis for the most gen- 
eral choice of complex phases. The predicted correlation 
between sin 2 #i 2 and sin 6*13 is shown in Fig|3] We 
find that, in the normal ordering case, sin#i 3 > 0.12 
is required to accommodate the non-maximal 1 — 2 
mixing. Besides the opportunity to measure this value 
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of #13 already in the Double Chooz experiment, this 
scenario is also promising for future searches of leptonic 
CP violation [2J. In the inverted ordering case, instead, 
there is no lower nor upper bound on sin #13. However, 
when e is small we find 0.05 < sin #13 < 0.12 (gray 
region in the lower panel of Fig{3|. We checked that 
this is the case when 771,3 0-05 eV, in analogy with the 
CP-conserving case (gray region in the lower panel of 
Fig§. 

Conclusions. We studied the most general mass matrix 
for two mass-degenerate neutrinos. This is possibly a 
good limit to understand the smallness of the 'solar' mass 
splitting Ato 2 2 relative to A77123. We have shown that 
such a mass degeneracy requires the equality of the 11 
and 22 matrix elements, as well as the vanishing of one 
CP violating phase. The first requirement points to a 
non-Abelian family symmetry, the second indicates that 
CP violation in the lepton sector should be not generic, 
if present at all. The matrix structure leading to the 
mass degeneracy is most easily accommodated, in the 
framework of three families, if neutrinos only mix in the 
1 — 2 sector and, therefore, the large 2 — 3 lepton mixing 
comes from the charged lepton sector. 

We have shown that all these features can be explained 
by the simplest quaternion family symmetry Q, together 
with the requirement of no CP violation. We discussed 
several realizations of our Q model, either employing sev- 
eral Higgs doublets or heavy flavon fields, and realizing 
the seesaw by either Higgs triplets or right-handed neu- 
trinos. In the limit where Am] 2 = 0, the model predicts 
also &i3 = 0. 

We discussed the possible perturbations generating 
nonzero Am] 2 and the consequent correlation with the 
nonzero value of #13. We studied in detail the predictive 
case where TOj, depends only on three real parameters 
plus one small perturbation. Both normal and inverted 
ordering of the mass spectrum can be realized. In the 
normal case, sin 6*13 > 0.12 should be found, close to 
the present upper bound. In the inverted case, sin #13 is 
smaller than about 0.1. A lower bound (sin #13 > 0.05) 
holds in the inverted case when the absolute neutrino 
mass scale 7773 (or equivalently the neutrinoless 2(3 effec- 
tive mass 777 ee ) is larger than about 0.05 eV. 

If CP violating phases are present, the Q model 
predicts a correlation between nonzero $13 and the 
deviation of #i 2 from the maximal value. Also in this 
scenario a lower bound sin #13 > 0.12 holds in the normal 
ordering case. 
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Appendix: Alignment of VEVs 

The Q model defined by eq.([6| contains three Higgs 
doublets and three Higgs triplets. Here we will discuss 
some features of the scalar potential responsible for gen- 
erating their VEVs. 

Let us consider the most general Q-invariant potential 
for (0i,02) ~ 2 (notice that the conjugate Q-doublet is 
(02, -#)): 

V = + 4<f> 2 ) 

+ §Ai[f> t 1 <M 2 + (0 2 2 ) 2 ] + A 3 (0Ui)(4<K) (20) 

+\H(4>[^) 2 + (4M 2 } + A 4 (0 t 1 2 )(40i) • 

It is bounded from below as long as Ai > | A 2 1 + | A3 + A 4 1 . 
For 777 2 < 0, the electroweak symmetry is broken in a 
minimum where |ui| 2 = |v 2 | 2 = — to 2 /(Ai + A3+A 4 — |A 2 |). 
The vacuum with v\ = and v\ — —ttt, 2 /Ai is a saddle. 
In order to make it the minimum, one needs to break 
softly Q to Z4, by adding 

V soft = m 2 (02>2 - <f>\<h) ■ (21) 

Notice that no other soft terms are allowed by the resid- 
ual Z4 symmetry. Then the electroweak symmetry is 
broken if 771 2 < |/x| 2 and the vacuum with v\ = and 
v\ = —(to 2 + /i 2 )/Ai is the absolute minimum as long as 
/^(Ai + A 3 + A 4 - |A 2 |) < to 2 (Ai - A 3 - A 4 + |A 2 |). This 
alignment of VEVs is the one analyzed in the text. 

In fact, the analysis would be pretty much the same, 
even if one does not want to introduce V so ft ■ I n this case 
the minimum is given by v± = u 2 = v (strictly speaking 
the relative phase between v\ and t> 2 can be ±1 or ±i, but 
for the sake of brevity we illustrate only one possibility) . 
Then, the charged lepton mass matrix in eq.([7]) can be 
rewritten as 

/ 71 75 ° \ / -vW \ 
m = \ -72 72 ^ v ^y* v ■ 

V 1 / V y 4 v 3 J 

Therefore, up to a maximal 1 — 2 rotation on the left and a 
harmless swap of ef and e 2 on the right, one recovers the 
same structure of to; as in the case V\ — 0. Notice that 
the maximal 1 — 2 rotation does not modify the structure 
of the neutrino mass matrix m u in eq.([7]), so that all 
predictions for the lepton mixing angles are unchanged. 
However, such rotation introduces FCNCs involving the 
electron, so that the non-standard Higgs bosons should 
be heavier to evade the experimental constraints. 

The addition of the third Higgs doublet ^3 ~ 1 does 
not change the qualitative features of the scalar potential 
discussed above. 
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Let us briefly discuss the origin of the VEVs of the 
three Higgs triplets in eq.(|6|, as well as of the two triplets 
(A 4 ,A 5 ) ~ 2, needed to generate m v in eq.([9]). The Q- 
invariant scalar potential reads 

Va = ELi MfA\At + Mf (A 4 A 4 + AiA 5 ) 

+ [/ilAi0i0 2 + H2A 2 {(j)l4>l - 4>2<i>2) + [i 3 A 3 (f> 3 (f) 3 (23) 

+ M4 (A 4 0i0 3 + A 5 020 3 ) + h.c] + Vl uarUc . 

In the hypothesis that Mi is much larger than the elec- 
troweak scale, one can integrate out Aj using its equation 
of motion. The triplet VEVs are thus determined as fol- 
lows: 

ul = -j^v lV2 , u* 2 = -^{vl-vl) , 



In the case v\ = 0, one has automatically w 4 = as 
assumed in eq. (|9j , however also u% vanishes, which is not 
acceptable. There are two ways to cure this problem. 
The first is once again to add terms which break Q softly, 
e.g. Vl° ft = Mf 2 (A\ A 2 + A 2 Ai), that allows to induce 
both u\ and u 2 from v 2 only. 

A second possibility without soft breaking is to resort 
to Q-invariant quartic couplings among Higgs doublets 
and triplets. Consider, in particular, 

X(4^ - 4<fc)(A| Ai + AlA 3 ) C Vl uarUc . (25) 

Then the equation of motion for Ai leads to u\ — 
Xu%\v 2 \ 2 / Ml . This contribution is tiny if all triplet 
masses are of the same order, but when Mi is instead 
close to the electroweak scale one can obtain u\ ~ u 3 . 
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